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Introduction 

A key ingredient in the design of engineering components and structures under general 
thermomechanical loading is the use of mathematical constitutive models (e.g. in finite 
element analysis) capable of accurate representation of short and long term stress/defor- 
mation responses. In addition to the ever-increasing complexity of recent viscoplastic 
models of this type, they often also require a large number of material constants to 
describe a host of (anticipated) physical phenomena and complicated deformation mech- 
anisms. In turn, the experimental characterization of these material parameters consti- 
tutes the major factor in the successful and effective utilization of any given constitutive 
model; i.e., the problem of constitutive parameter estimation from experimental mea- 
surements. 

Traditionally, simple, basically trial-and-error procedures (graphical/mechanistic fitting) 
have been used for simple models, but these are certainly rather limited in more general 
situations Fig. 4 . This is particularly true in dealing with very large number of material 
constants that are often lacking in their direct physical interpretation, where complica- 
tions due to the vastly different scaling and highly interactive nature of these parameters 
in a large test matrix under various controls (stress, strain, or mixed) under transient and 
steady-state conditions Fig. 2. 

An urgent need and obvious need therefore exists for a systematic development of a 
general methodology for constitutive parameter estimation. This provides the main moti- 
vation for the present work. 

Background and Approach 

The problem belongs to the class of inverse problems [1] of mathematical programming 
and optimization theories. Its solution requires three major and interrelated parts Fig. 4 
in its application to the present dynamic (time-variant) case; i.e., (a) primal analysis 
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(response functionals) for the differential form of the constitutive model, (b) sensitivity 
analysis, and (c) optimization of an error/cost function. The optimization algorithms for 
the last part (c) are presently very well developed [2]. This is not the case, however, 
regarding work on the other two, intimately related, parts (a) and (b). Mainly due to the 
greater mathematical complexity and associated intensive computational demands for 
the present dynamic and nonlinear case, compared to other more traditional optimization 
problems of linear structures. This renders unsuitable or even inapplicable several of the 
available solution methods and algorithms for primal and sensitivity analyses, Fig. 5. 

For example, using an explicit integration method, with its known material-dependent 
conditional stability limits, becomes very ineffective for primal analysis, in which essen- 
tially thousands of “different materials” are being processed for response predictions (in 
a typical optimization cycle), thus making any adaptive time-stepping strategy very com- 
plicated, if at all possible. Similarly, in addition to several accuracy and numerical insta- 
bility problems, the use of finite-differencing schemes for sensitivity analysis can easily 
become computationally prohibitive with the increase in material constants and time win- 
dows for fitting with long-duration tests. 

The highlights of the mathematical formulations and main features of the present devel- 
opment (Figs. 6 and 7) are summarized as follows. Posed as a least-square, con- 
strained, nonlinear mathematical optimization [2], we use an objective functional of the 
minimum-deviation-error type, i.e., differences in the predicted and measured responses 
at varying times. The material constants constitute the design variables, with several 
(side) constraints to ensure a physically-meaningful model. For the primal response 
analysis part, we utilize an implicit, unconditionally-stable, integration algorithm. Details 
and several applications of this scheme are described in a separate presentation in this 
proceeding; see also [3]. The sensitivity and analysis is of the direct type preformed on 
the basis of an explicit, recursive, form associated with the above integrator. Finally, the 
optimizer segment is of the gradient-based type, utilizing a sequential-quadratic pro- 
gramming scheme [4]. It is these three approaches combined that provide for the robust- 
ness and computational efficiency. 

The overall strategy is summarized in the flow chart of Fig. 8. Its main driver (dubbed 
COMPARE for Constitutive Model B^Bamter Estimator) controls the three solution mod- 
ules (primal analysis, sensitivity, optimizer), together with the management of data files 
and results. From the practical standpoint (Fig. 9), the overall strategy is sufficiently gen- 
eral to handle comprehensive test-matrix data, under arbitrary load-control variables, 
multiaxial stress/strain tests, and transient as well as steady-state response measure- 
ments. 

Results 

All applications given utilize a viscoplastic model of the nonlinear kinematic-hardening 
type, GVIPS [5], having a total of eight material constants (Fig. 10). The first part of the 
studies [Figs. 11-15] are performed on a simulated (“exact") model material, for valida- 
tion and to investigate the issues of parameter sensitivities, accuracy, comparisons, and 
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In the second part [Figs. 17-22 J # we use actual test-matrix results for TIMETAL21 S 
material. This includes three tensile test under different strain rates, creep tests with 
three different imposed stresses, a relaxation test, as well as a three-step creep test. 
The “fitting” success in this latter, more realistic application, with vastly different condi- 
tions and very large number of data points, clearly points to the potential benefit and 
practicality of the general methodology. 
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MOTIVATION 

Two major obstacles to fully utilizing recent time-dependent/hereditary con- 
stitutive models in practical engineering analysis: 

• Lack of efficient and robust integration algorithms 

• Difficulties associated with characterizing large number of required 
material parameters 

- Most material parameters lack obvious/direct physical interpretations 

- Even under load histories in simple laboratory tests, several parameters will 
highly interact to affect predicted responses 

- Further complications due to: 

(i) Incompleteness of response measurment in both time and state 

(ii) Vastly different scaling of constitutive parameters 

- Urgent need exists for specific guidelines in the systematic development of 
general methodology for constitutive- parameter-identification 


Fig. 2 
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OBJECTIVES 

General: Systematic development of general methodology for constitutive 
parameter estimation 


Specific: 

• Mathematical formulation of a basic optimal material parameter estimation 
scheme 

• Computational algorithms for implementation 

• Validation tests and performance studies 

• Alternatives for further refinements in fitting; 

that is, variable/cost function scaling, weights, multicriteria optimization theo- 
ries, etc. 


Fig. 3 


BACKGROUND 

• ‘Traditional” approaches for constitutive parameter estimations: 

- Essentially trial -and-error in nature 

- Based on several assumptions about test conditions and material behavior 
that are rarely satisfied in actual tests 

- Difficult to control error propagation in sequential evaluations of parameters 

- Rather limited in applications 

• Modern approaches for constiutive parameter estimations based on 
mathematical programming and optimization theories: 

- More aeneral/ svstematicallv derived 

- Three major parts for modular implementations: 

(a) Primal analysis for response functionals 
(integrated history) 

(b) Sensitivity analysis 

(c) Optimization of “cost” function 

Fig. 4 
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BACKGROUND 


* Methods and algorithmic details differ greatly for the sensitivity analysis: 

- Finite difference methods (prohibitively expensive; prone to errors) 

- Evolutionary sensitivities approach (expensive two-subproblem 
integrations; special coding for number of parameters dependent sizes of 
arrays; stiffness/singularity problems) 

- Adjoint sensitivity approach (regressive computations with large storage 
requirements; or increased computational cost for terminal adjoint 
problems) 

- Direct-differentiation sensitivity (most effective and accurate when 
consistently derived with the underlying implicit integration of the model) 
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MATHEMATICAL FORMULATION 


• Framework required characteristics: 

- Coupled nonlinear system (internal/external state variables) 

- T ransient response with different possible steady-state conditions 
(time variance) 

- Arbitrary (optional) control variables; i.e., stress-, strain-, mixed-types of 
loading (general test matrix) 

• Approach: 

- A least-square, constrained, nonlinear mathematical optimization problem 

- Material parameters constitute the design variables, with several associated 
side constraints to ensure physically meaningful! model 

- The technique is of the minimum-deviation-error type, for the integrated mul- 
tiaxial response (functionals) 


Fig. 6 
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MATHEMATICAL FORMULATION 


• Noteworthy aspects: 

(i) Robustness and effectiveness: 

- Unconditionally stable implicit integration for primal analysis (model 
problem) 

- Direct -differentiation approach for accurate sensitivity analysis 

• State-of-the-art optimizer using sequential quadratic programming (with 
exact gradient and variable metric/Hessian) for least-squares 
minimization 

(ii) Computational efficiency: 

- Non -iterative “ exact " sensitivities (once after primal analysis time-step 

- Effective scaling for both the design variables and the objective function 
in optimization 

- Enhanced iterations with line searches in implicit integration 

*9.7 


COMPUTATIONAL ALGORITHM 



Sequential Quadratic 
Programming (SQP) 



Direct Differentiation 
Approach 


Implicit Integration 
For Primal Analysis 


/(COMPARE)^ 
V Driver J 


- Identify active/passive variables 
for a test 


- Scale design variables and 
objective function 

- Formulate a single design 
optimization problem 

- weighted objective function 

- constraints 

- sensitivities 


Results 


Fig. 8 



Analysis Data: 

- Problem Type/Control 

- Muitiaxial Responses 

- Time Window 

Estimator Data: 

* Number of Tests 

- Initial Design Variables 

- Upper/Lower Limits 

- Active/Passive 
Variables 

- Variable Grouping 

- Weighting Factors 

Optimizer Data: 

- Convergence 
Tolerance 

- Iteration Limits 

- Stop Criteria 
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COMPUTATIONAL ALGORITHMS 


Features and capabilities: 

- General test matrix with arbitrary control (stress, strain, or mixed) with 
multiaxial measurements (two normal plus one shear component), e.g., 

Stress-control Creep 

Mixed-control Relaxation 

Stress-control Tension/constant stress rate 

Mixed-control Tension/torsion test 

- Active/passive design parameter activation (e.g., parametric study) 

- Upper/lower side constraints 


- Arbitrary number of tests and time windows selected for fitting 

- Goodness-of-fit statistical measures 


Fig. 9 


SENSITIVITY ANALYSIS 

* Model problem 

- Unified viscoplasticitv with potential GVIPS (nonlinear hardening: static 
recovery mechanism; isotropic/nonisotropic yielding) 

- Total of seven viscoplastic parameters (design variables) 

Flow equation = 2 (viscosity p; exponent n) 

Evolution equation =4 2 (hardening: modulus H, exponent p) 

+ 2 (recovery; modulus R, exponent m) 

Yield threshold = 1 (k^) 

• Simulated “actual” material (perfect model representation capability) 

- Normalized sensitivity plots 

- “Accuracy” and efficiency comparisons with traditional finite-differencing 
schemes 

- Creep, relaxation, constant-strain-rate tension as representative tests 

Fig. 10 
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SENSITIVITY ANALYSIS (SIMULATED TESTS) 


DESIGN 

PARAMETERS 

T, 

T* 

T 50 

IMPLICIT 

F.D. 

IMPLICIT 

F.D. 

IMPLICIT 

F.D. 

K 

1.7899 

(10)-1 

1.7864 

(10)-1 

4.5844 

4.5816 

4.5844 

4.5817 

n 

-1.1361 

(10)-1 

-1.1384 

(10H 

2.1588 

(10)-1 

2.1525 

(10)-1 

||g| 


n 

1.4583 

(10)-6 

1.4570 

(10)-6 

-1.7294 

<10)-6 

-1.7291 

(10)-6 



m 

9.2111 

(10)-19 

0.000 

-6.3965 

(10)-6 

-6.4033 

(10)-6 

-1.5127 

(10)-5 

-1.5144 

(10)-5 

p 

4.2054 

(10)-3 

4.2041 
(1 0)-3 

-9.4439 

(10)-1 

-9.4236 

(10)-1 


| 

R 

-91043 
( 1 0)-9 

0.000 





H 

1.4634 

(10)-7 

1.4628 

(10)-7 

3.6125 

(10)-5 

3.6085 

(10)-5 

3.6124 

(10)-5 

3.6084 

(10)-5 


► Time Evolution 


Fig. 11 


SENSITIVITY ANALYSIS (SIMULATED TESTS) 

TIME EVOLUTION 

16 x 10 s 
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Relative (least-squares) Error 


EFFICIENCY OF IMPLICIT SENSITIVITY SCHEME 



TENSILE 

CREEP 

RELAXATION 

RELATIVE 

CPU 

IMPLICIT 

1.0 

1.0 

1.0 

F.D. 

8.004 

8.235 

7.931 


General Estimate: 

Relative CPU a&xN.xN^ xN.xt\ 

l TV l 

Nj = Number of integration steps per each (equal) optimization 
interval 

N w = Number of (equal) sampling points (time windows) per test 
N t = Number of total tests 

t) = Trial-Differencing-Accuracy Factor q = 2 10 -» ? 

Fig. 13 


TYPICAL CONVERGENCE AND FITTING ACCURACY 

(SIMULATED MATERIAL FOR TENSILE/CREEP/RELAXATION TEST) 



Number Of Iterations Time (sec) 


Convergence Accuracy in Tensile Test 

Fig. 14 
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TYPICAL CONVERGENCE AND FITTING ACCURACY 

(SIMULATED MATERIAL FOR TENSILE/CREEP/RELAXATION TEST 

180 r 



Accuracy in Creep Test 

Fig. 15 


160 - 
140 - 
120 - 
100 - 
80 - 
60 * 

40 ■■ 

i 

20 

0 9000 18000 27000 36000 45000 

Time (sec) 

Accuracy in Relaxation Test 


x True Set 

— Initial Set 

— Optimized Set 


APPLICATIONS 

• Material: -TIMETAL21S 

- Temperature = 650°C 

• Experimental tests available 

(a) 3 tensile tests 

(b) 3 creep test 

(c) 1 relaxation test 

(d) single 3-step creep test 

• Results/Studies 

- varied number of tests included in fitting 

- varied sampling-time intervals within each test 

- varied material-parameter bounds and optimization weights 

- comprehensive case: all tests (a)-(d) included in fitting 


Fig. 16 
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SUMMARY OF COMPREHENSIVE FIT CASE 

( TIMETAL 21 S; 8 TESTS; T=Tensile, C=Creep, R=Relaxation, SC=3-Step Creep) 



T1 

T2 

T3 

Cl 

C2 

C3 

R 

SC 

Number of 
Fitting Points 

93 

67 

57 

96 

103 

107 

72 

285 

Weight Factors 
(equal weight) 

0.125 

) 

0.125 

0.125 

0.125 

0.125 

0.125 

0.125 

0.125 

Weight Factors 
(variable weight) 

0.1193 

0.1657 

0.1947 

0.116 

0.1077 

0.1037 

0.1542 

0.0389 


a) Fitted Points and Weights 



Number of 
Iterations 

Number of 
Function Calls 

Number of 
Gradient Calls 

Normalized 

CPU 

Estimated 

(equal) 

53 

61 

54 

1.0 

Estimated 

(variable) 

49 

54 

50 

0.882 


b) Solution Efficiency 

Fig. 17 


SUMMARY OF COMPREHENSIVE FIT CASE 

{TIMETAL SIS, 8 TESTS) 



Fig. 18 
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TENSILE BEHAVIOR CORRELATION 

(COMPREHENSIVE FIT) 

Given: 

- a “far" initial guess 

- with/without bounds on parameters 
unique “optimal” response 

,-6 


6q = 8.333x10 


x Experiment 

Initial Set 

Optimum Set 

x (Rcieased-Constraints) 



400 800 

Time (sec) 


1200 




0 20 40 60 80 100 120 

Time (sec) 


2 4 6 8 10 12 

lime (sec) 


Fig. 19 


CREEP BEHAVIOR CORRELATION 

(COMPREHENSIVE FIT) 


Given: 

- a “far” initial guess 

- varying sample-time intervals 
unique “optimal” response 



0.04 
0.03 1- 


$ 0.02 b 


CT 0 = 72.4MPa 



3000 6000 9000 

Time (sec) 
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or 2 2(MPa) 


RELAXATION BEHAVIOR CORRELATION 

(COMPREHENSIVE vs. SINGLE-RELAXATION TEST FIT 


Given: 

- a “far” initial guess 

- variable optimization weights 
-> unique optimized response 



Time (sec) 

(a) Comprehensive Fit 


Comprehensive fit gives better overall 
predictions with comparable accuracy 
to the “pure” single-response curve 
correlation 



Time (sec) 

(b) Single-Test Fit 

Fig. 21 


3-STEP CREEP CORRELATION 

(COMPREHENSIVE FIT) 



0 1 2 3 4X10 4 

Time (sec) 

Fig. 22 
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SUMMARY/CONCLUSIONS 


• Validation with exact (simulated) material 

- Given an accurate constitutive model, exact correlation is achievable by 
COMPARE 

• Assessment with real materials (TIMETAL 21 S) 

- Automated material parameter estimation enables the model to achieve its 
“ best " correlation 

- Serves as a tool for identifying critical experiments to maximize pertinent 
“data content” (e.g., one test for tension, creep, cyclic and (initial) relaxation) 


- Requires minimum “user expertise” 

- Gives a measure for model suitability and directions for its further 
enhancement in realistic range of applications 

- Estimates for effects of model versus experimental (noisy data) deficiencies 
(COMPARE'S knowledge of | 


Fig. 23 


SUMMARY/CONCLUSIONS 


- Including more data points in fitting enhances the optimizer convergence 
speed 

• New model implementation 

- Demands mfiCfe than just a definition of elementary (differential) flow/evolution 
equations - integrated form Jacobian 

Fig- 24 


FUTURE WORK 

• Inclusion of new material models in COMPARE library 

• Experience in applications with multiaxial test fitting 

• "User-friendly” enhancements 


Fig. 25 
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